Spin-disordered superfluid state for spin-1 bosons with fractional spin and statistics 
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We study a strongly correlated spin-1 Bose gas in 2D space by using the projective construction. 
A spin-disordered superfluid state is constructed and proposed as a candidate competing with the 
conventional polar condensate when interaction is antiferromagnetic. This novel state has a non- 
trivial topological order whose low energy excitations carry fractional spin, charge, and statistics. 
The spin excitations become gapless only at the edge and are described by level-1 SU(2)spin Kac- 
Moody algebra. The edge state is identical to the edge state of the chiral spin liquid or the right 
moving sector of spin-1/2 chain. 
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So far most studies of ultra cold alkali atomic gases 
have been focused on Bose-Einstein condensation (BEG) 
in weakly interacting dilute limit (for review, see for in- 
stance [Q). By contrast, a strongly correlated Bose gas 
may lead us to interesting novel phenomena. While it 
is unclear whether one can achieve a stable strongly cor- 
related gas in 3D by simply cranking up the scattering 
length without collapsing the system, a 2D gas with even 
relatively weak interaction can be strongly correlated in 
nature. The argument is based on the renormalization 
group analysis that showed that the interaction in 2D 
is marginally irrelevant only in a dilute limit specified by 
lnln(l/na^) 3> 1, where n is the particle density and a 
can be thought of the interaction range or the scattering 
length. The double logarithm imposes a much stronger 
condition of the validity of the dilute limit in 2D than the 
familiar na^ ^ 1 in 3D. For generic BEG system such as 
the disk condensate of ^■^Na of Ref. Hi or a microelec- 
tronic chip of condensed ^^Rb atoms Ml, lnln(l/na^) is 
of 0(1). These systems are, at least, not weakly corre- 
lated. Now the question is that what is the ground state 
for those not-weakly correlated 2D boson gas? In this 
letter, we shall show that a new class of 2D superfluid 
can emerge due to the strong correlation. 

a. Proposed spin- disordered superfluid. Let us con- 
sider a strongly correlated gas of spin-1 bosons (pm {m = 
0, ±1) in a 2D homogeneous space with generic two-body 
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J^'^* with 
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where C^fm^ = (S'5'z|ll; mim2) are the Glebsch-Gordan 
coefficient for total spin S = 1 + 1. In the alkali 
atomic gases, the interactions for all three spin chan- 
nels Vs, S = 0,1,2 are commonly short-ranged. The 
Hamiltonian has a phase and spin rotation symmetry: 
U(l)chargcX SU(2)spi,i. The Hamiltonian derived by Ho 
Q, ^'"t = ^(con^ + CaS^) where n and S are the den- 
sity and spin operators, corresponds to a special form of 
(I) with Voir - r') = (co - 2c2)S{r - r'), Viir - r') = 
(co-C2)<5(r-r'), andy2(r-r') = (cq C2)<5(r - r'). The 
case of antiferromagnetic spin interaction (V2 > Vq or 
correspondingly C2 > 0) is somewhat intriguing. Mean 



field theories g realize that the (coherent or frag- 
mented) polar condensate is energetically favored for this 
case. We note that the polar condensate, no matter co- 
herent or fragmented, is a spin ordered state, as it directly 
condenses single bosons {{(pm) 7^ 0) and therefore always 
breaks SU(2)spin symmetry in the thermodynamical limit. 
For spin-1 bosons in an optical lattice, Demler and Zhou 
speculated possible spin disordered (liquid) phases and 
fractionalized spin excitations for arbitrary dimensions 
But, a question remains: Is there a spin-disordered 
superfluid competing with the polar condensate? 

We have found such a spin-disordered state. It can be 
visulized as a gas of closed (connected) clusters consti- 
tuted by various number of bosons; the size I can take up 
the value of 2 (dimer), 3 (triangle), . . ., up to iV for an 
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N boson system (see Fig. 0). Each cluster is a spin sin- 
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FIG. 1: A cluster configuration for = 11 spin-1 bosons with 
arrows indicating spin. 'Dotted line'= g{ri — Vj). 

glet and is totally symmetric among all bosons within 
the cluster. Mathematically, an Z-cluster is described 
by G({r„mJ,=i..,) = IlU - r.+i) tr{nL, e-^"'} 
(with r/_|_i = ri identified), where g(r) is a pairing wave- 
function of Pa;-|-ij,-symmetry (thus odd in r) , e is a 2 x 2 
antisymmetric matrix, and "^"^ is a (symmetric) 2x2 ma- 
trix defined by the Glebsch-Gordan coefficients for ^ ® ^ , 
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{{SSz)lm\^^; aP) . The state is a superposition of 
all possible cluster configurations specified by {ni} with 
ni the number of Z-clusters. The wavefunction for the 
component of the state with N bosons is, apart from an 
overall factor, 



xGJK,m?}ti), 



(2) 



where Nc — X); the total number of clusters, p is 
a label used to sort all Nc clusters, and Ip is the size of 
the p-th cluster. In this labelling, the boson spins and 
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coordinates {ri,mi}^^ are ono-to-one mapped onto the 
set {r^,TO^;p — l-.N^;! — l..lp}. The wavefunction is 
symmetrized for N bosons, as indicated by S^. 

Table. | compares the interaction energies of the two 
states. One sees that the spin-disordered state we pro- 
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TABLE L Interaction energies calculated at mean field level. 
Notations: no is the density of condensed bosons; Sr = r — r' ; 
Vs(5r)'s are defined in Eq. (^). Note that the last two inter- 
action energies for the spin-disordered vanish for 5r 0. In 
our calculation, the fermion representation Q is normalized 
by requiring {(j>ln,<l>m) equal to the boson density. 



posed has a fower interaction energy when V2 is positively 
large and short ranged. Another important feature is 
that wherever two bosons in the same spin state (except 
m = 0), say located at r and r', approach each other, ^'^^r 
vanishes as ~ (r — r')^, as implied in Table |. Readers 
who are not interested in the technical details may skip 
now to the end of the paper for other interesting physical 
properties. 

b. Projective construction: a theoretical technique. 
Unlike the usually studied, weakly interacting dilute 
limit, a Bose gas of strong repulsion cannot be treated 
by conventional perturbation theory. How to derive the 
groundstatc wavefuntion is obviously challenging. Fur- 
thermore, even if one has such a poposed wavefunction 
as the above, it is still difficult to extract the low 
energy excitation properties, since it is a highly 'entan- 
gled' state. The analogue of Gross-Pitaevskii equation 
does not exist yet. To make some progress towards such 
a difficult task, we borrow the 'projective construction' 
method from quantum Hall studies, which is a standard, 
successful approach in parallel with the Laughlin wave- 
function approach (see Ref. |^ and references therein). 
In this spirit, we introduce two spin-i fermions, ^aa, 
each with two (physical) spins a =t, i and two fictitious 
"colors" labelled by a, 6. The color is necessary to furnish 
a minimal spin-I bosonic representation at every space- 
time point. The boson can then be represented by 



0m (r) = ipaa{r)ipbi3{r)eab'<fap ■ 



(3) 



Note that the color degree of freedom is non-physical. 
All physical states or operators are required invari- 
ant under a local SU(2)coior transformation: ipaaii^) 
W{r)ab'4'ba{T^) , where W{r) is an SU(2)coior matrix. One 
can quickly check that the boson operator (jjm defined 
in (^ is indeed a color singlet, perfectly invariant un- 
der above transformation. Eq. ^ allows us to construct 
physical boson many-body wave function from the 
unphysical fermion many-body wave function j^'-'') by 
projecting it into the "color" singlet sector. In a mathe- 



matical equation, that means 

*^ (n , mi ; r2 , m2 ;...) = (0| n. 0m. (r. ) I , (4) 

whose explicit form is precisely ^'5^r of Eq. (^. 

The relationship (Q) suggests that the low energy ef- 
fective theory of our spin-1 boson system can either be 
described in terms of the boson operator or equiv- 
alently, in terms of the fermion operator ^aa- Unfor- 
tunately, for a strongly correlated 2D system, there is 
no known rigorous way to derive the effective theory 
from a microscopic Hamiltonian like (^. One usually 
first writes down a most natural form of it on sym- 
metry grounds, checks its stability against interactions 
in low energy limit, and finally compares it with ex- 
periments. In this spirit, the effective theory for the 
state ^''(ri, mi; r2, m2; ...) is described, in the fermion 
description, by a theory of independent fermions cou- 
pled to color SU(2) gauge fields. The gauge field is de- 
noted as Af^ = ^T^aJ^, I = 1,2,3, where are the Pauli 
matrices generating the SU(2)coior algebra. The gauge 
fields are introduced to project out the unphysical col- 
ored excitations.!^ The effective theory is then 

= 'i'lpLiDo)ab1pba + ^IpLiD ■ D)abi'ba 

+all symmetry allowed interactions , (5) 

where (D^)ab = 5abdf_, - i{Af_,)ab (for notation, see [^) 
are the covariant derivatives. From the fermion effective 
theory, we can study various fermion states, which, after 
the projection (|^), lead to various physical boson states. 

To see which fermion states are likely to appear as 
the ground state, we need to consider the interactions 
between the fermions. Interactions can be either origi- 
nated from the boson-boson interactions or dynamically 
generated by gauge interactions. As an example of non- 
Abelian gauge theory, the study of QCD |11 shows that 
the Yang-Mills gauge fluctuations can generate a strong 
attractive interaction between quarks due to the instan- 
ton effect, which leads to quark confinement. In our case, 
the gauge fields are used to mediate a strong attractive 
interaction between color-opposite fermions, since by def- 
inition (^) the gauge interaction is supposed to bind two 
fermions locally into a colorless boson. Therefore, the 
strong SU(2) gauge interaction naturally leads to color- 
singlet Cooper pairing. 

Let us consider two simplest color-singlet parings: 



(V'aa(r)V6/3(0)) 

i^.(|r|), 



(s, spin triplet) 



\ Eafc Cap Rp{\r\) {x + ijj) , {px+iy, Spin singlet) 

where m can be or 1. Both s- and Px+iy-'^&ve states 
produce a full gap on the Fermi surface. i?s,p(|r|) are 
complex functions of |r|, generically expected to mono- 
tonically fall off exponentially at large distance. 

c. s-wave pairing: conventional BEC phases. Con- 
ventional BEC phases are easily recovered through the 
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s-wave paring channel. In the limit of strong confine- 
ment, i?s(|r|) becomes a delta function ~ S{r). 

Polar condensate. This is a special kind of spin 
nematic state. Fermion confinement occurs in spin-1, 
m — channel. Here, the pairing in Eq.(^) reduces to 
{tpaa{ic)ipbi3{ic)) = ^/pe"^ '^ab'^^f'^ where one may think 
of p related to the condensed boson density. The result- 
ing state is nothing but the so-called polar condensate 
with ((/i,„=o) 7^ 0. This state breaks both U(l)chargeaiid 
SU(2)spin invariance, and was considered to be favored if 
the spin interaction is antiferromagnetic (V2 ^ Vq) 

Ferromagnetic condensate. This case is the same 
as the polar condensate except that the fermions are 
confined into the m = 1 (or equivalently —1) chan- 
nel. The order parameter becomes (V'aa(r)'0()/3(r)) = 
\fpe^'^ ^ah^™fi^^ which corresponds to the ferromagnetic 
condensate. Like the polar condensate, it breaks both 
U(l)chargoaiid SU(2)spijj invariance. This state is pre- 
sumably favored if the spin interaction is ferromagnetic, 
1^2 < 14). 

d. Px+iy -state: topological superfluid. The spin- 
disordered Px+iy-vfave pairing energetically competes 
with the polar condensate for a large positive V2 (see 
Table. ||). Assuming the state to exist, its low energy ef- 
fective theory can be routinely constructed, simply based 
on the broken (physical) symmetries without relying on 
microscopic details. We find the effective Lagrangian of 
the state. 



ka(r)A(r, r')e./5[e^/^' ^■'''tT\,i;l{r) + h.c 



(7) 



where A(r, r') is the p-wave pairing wavefunction. In the 



momentum space it becomes A(k) 



'-{kx + iky). The 



theory is evidently invariant under the global SU(2)spi„ 
and local SU(2)coior transformations. However, since 
the Px+iy-wa.ve condensate (see Eq. (^) spontaneously 
breaks the charge U(l), the phase fluctuations of A gives 
rise to the usual (gapless) superfluid mode. The charge 
sector therefore is rather conventional. In contrast, we 
shall show that the spin sector has excitations with frac- 
tional spin and statistics. 

Let us focus on the spin sector of the effective the- 
ory (^, by regarding A as a fixed non-dynamical field. 
To see the symmetry more clearly, we introduce four 2- 
component spinors, each uniting a pair of particle and 
hole operators: 

riac^f t V (a,6=l,2;a,/?=t,i). (8) 

In terms of 77's, the effective action in momentum space 
becomes 



&{k) 
hik) 



= -i / d^k4jk)Mk^^ - A^')U;b0Vbpik) , (9) 

= f° ® r° ® [^iv + h{k) • a] , A:^ = {lj, k) , (10) 
= {Aoky, Aokx, kV2Af-/^) . (11) 



The !^{k) is better viewed as an 8 x 8 matrix in an or- 
thogonal basis {f" = l,f'}spin® {r° = 1, Tricolor <8){cr° = 

1,(7'} {I = 1^,3) with the cr's acting on the "Nambu" 
spinor space (||). A careful reader may recall = ^t'o^. 

The spectrum of fermionic excitations is determined 
by the pole of the inverse of S'{k) at = 0, E(k) = 



(k2/2M - fi) + Al . The fermions are fully gapped, 
so we integrate out them to find the effective action of 
gauge fields. To do so, we have chosen a background field 
gauge in (^ such that = Const but non-commuting, 
^ 0. The approach we adopted here is stan- 
dard in non-Abclian gauge field theory. Since the action 
(^) is quadratic in fermion fields, the gauge effective La- 
grangian is exactly given by 

■^cs[A] = -^^Jd^kTi\n^{k-A) (12) 

where the 'Tr' is over the internal space of f®T®(T. With 
the matrix given in ( p^ ) , a straightforward calculation 
gives 



■^es[A] = %^e^^P\x{A^d,Ap - fA^A.Ap} 



(13) 



where Q is a (topological) winding number {Q — I [p^) 
and the stands for higher derivatives, including a 
Maxwell term. Now the 'tr' in ( p^ is over only the color 
gauge indices. The effective theory of gauge bosons is 
thus a level 1 (non-Abelian) SU(2) Chern-Simons the- 
ory. All gauge excitations gain a dynamically generated 
topological mass (i.e., gapped). The gauge interaction 
becomes short ranged. So the effective theory we found 
in (0) is stable. 

e. Physical properties of the Px+iy-state. After ob- 
taining the low energy effective theory, we are ready to 
study the measurable properties of the state. First the 
Px+iy-state is a superfluid which does not break spin ro- 
tation symmetry. The only gapless excitation is the su- 
perfluid mode described by the phase of A. All spin 
excitations have a finite energy gap. Due to the Chern- 
Simons term, the SU{2) gauge field is not confining. 
Thus the excitations described by ■0 (or 77) have a fi- 
nite energy gap (instead of infinite energy gap). Those 
excitations carry spin-1/2 and one half of boson charge! 
They also have a semion statistics (i.e., the statistical 
angle is 9 = 7r/2, right between boson and fermion), 
as implied by the level-1 SU{2) Chern-Simons terms. 
The spin rotation symmetry implies that {4>m) = 
for all TO. However, two-boson and three-boson opera- 
tors both can have finite expectation values. In short 
distance, {(j>m(zi)4im'{z2)) ~ -4\/3(zi - Z2)'^C'^^, and 

{(f>rniizi)(l)m2iz2)4'm3{z3)) ^ 8V2(zi - Z2)(z2 - ^3)(23 " 



Zl)er. 



where z = x-\-iy. The spin-disordered super- 



fiuid has an unusual off-diagonal long range order. We 
see that the minimum vortex has one unit of quantized 
vorticity. Our Px+iy spin-disordered state breaks the par- 
ity and time reversal symmetry. The total angular mo- 
mentum of the ground state is h per boson. (Note such 
a total angular momentum is equal to the total angular 
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momentum of usual boson superfluid with one vortex at 
its center.) Thus, spinning the bosons may help to create 
the Px+iy spin-disordered superfluid. 

/. Edge excitations In the rj bases, the mean-field 
fermion Hamiltonian (described by (H) with being 
set to zero) contains four identical 2x2 blocks: H = 
/i(k) • a. In each block, the function ft,(k) defines a 
mapping from the k-space to with a winding num- 
ber 1. This non-trivial winding number leads to a unit 
Hall conductance. JT^, Q Since each block contributes 
a unit Hall conductance, it leads to an edge state sim- 
ilar to the one from v — 1 quantum Hall state, as 
required by gauge invariance.|]l6|, p^ , 18 Such an edge 
state can be described by one chiral fermion field Xaa 
(one for each block). Therefore, if we ignore the SU{2) 
gauge fluctuations, the mean-field gapless edge excita- 
tions of our px+iy-state are described by the following 
effective theory: Xl^^i{dt — vdx)Xaa- Only two of the four 
A's are independent because rjaa (see (ph) are Majorana 
fermions and each Xaa is obtained as a linear combina- 
tion of the two components of spinor rjaa- Therefore, 
the mean-field edge state contains only two independent 
branches of chiral fermions. Obviously the mean-field 
edge effective theory has SU(2)coior and SU(2)spin sym- 



metries generated by t 



^0 



and 



I = 1,2,3, 



respectively. Having the SU(2)coior x SU(2)spin symme- 
try and a central charge c = 2 (i.e., two branches of 
chiral fermions), we find that the mean- field edge state is 
described by SU(2)coior xSU(2)spin Kac- Moody current 
algebra of level 1. After including the SU(2)coior gauge 
fluctuations to go beyond the mean-field theory, the edge 
effective theory becomes 



Xi^i m - iAo) - v{dx - lAx)] 



ah ' 



(14) 



The effect of SU(2)coior gauge fields is to remove the 
SU(2)coior sector of the Kac-Moody algebra from the 
low energy spectrum. ||l^ Thus the physical edge state of 
the px+iy-siaie is described by a level-1 SU(2)spin Kac- 
Moody algebra. The physical edge state is identical to 
the right moving sector of spin-1/2 chain. Despite the 
finite gap in the bulk, the spin excitation is gapless at 
the edge. The operator that creates the gapless spin-1/2 



quasiparticle on the edge is given by the spin-1/2 pri- 
mary field Va{x,t) in the SU(2)spin Kac-Moody algebra 
which has a scaling dimension 1/4. The quasiparticle 
propagator has a form (ya{x^t)Va{Q)) ~ (x — vt)~^/'^. 
The boson operator 0„i on the edge becomes the spin-1 
primary field which is the spin current operator on the 
edge. The boson propagator on the edge is given by 
{<j)m{x, t)(f)l^{Q)) '-^ [x — vt)~'^. (The boson propagator is 
short ranged in the bulk due to the finite spin gap.) This 
will lead to a non-linear I-V curve / oc for boson 

tunneling between two edges. The spin-1/2 quasiparticles 
can tunnel between two edges separated by a bulk Px+iy- 
state. The tunneling I-V curve has a form / (x 
in the weak tunneling limits. Finally, we briefiy mention 
that the polar condensate has gapless spin excitations in 
the bulk whereas the Px+iy spin-disordered superfluid has 
a gapless spin excitation only at the edge. A dramatic 
difference can be seen in the spin susceptibility by NMR 
experiments. 

g. Conclusions A two-dimensional boson gas in 
ultra-cold alkali atomic systems can be strongly corre- 
lated. A spin-1 boson system can have a very inter- 
esting spin-disordered superfluid state, which carries a 
non-trivial topological order. |]l9) Such possibility is inter- 
esting, since they might exhibit some of the spin liquid 
phases that have been long theoretically speculated in the 
context of high Tc superconductors but never been clearly 
identified by experiments. We believe that the alkali 
atomic gases may provide the first controlled laboratory 
to check those speculated theories and enrich our under- 
standing of the strongly correlated systems. In fact the 
Pa:+ij^-state is closely related to the chiral spin state. pO[| 
The two states have the same bulk effective theory de- 
scribed by level-1 SU(2)coior Chern-Simons theory, |2l| 
and the same edge effective theory described by level- 
1 SU(2)spin Kac-Moody algebra. |16[ The spin sector of 
the two states are identical. 
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